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Contrary to expectations, several cosmic ray events with energies above 10^° eV 
have been observed. The flux of such events is weU above the predicted Greisen- 
Zatsepin-Kuzmin cutoff due to the pion production (via the A resonace) of ex- 
tragalactic cosmic ray protons off the cosmic microwave background. In addition 
to the relatively high flux of events, the isotropic distribution of arrival directions 
and an indication of hadronic primaries strongly challenge all models proposed 
to resolve this puzzle. Models based on astrophysical accelerators need to invoke 
strong Galactic and extragalactic magnetic fields with specific properties which are 
yet to be observed, while models based on physics beyond the standard model of 
particle physics generally predict photon primaries contrary to experimental indi- 
cations. The resolution of this puzzle awaits a significant increase in the data at 
these energies which will be provided by future experiments. 



1 Introduction 

The role of cosmic rays in the discovery of fundamental properties of nature 
dates back to the early 1900's. Over the last century, Brazilian physicists have 
actively contributed to these discoveries with the most notable example of Prof. 
Cesar Lattes and the discovery of the tt meson decay in Chacaltayald After a 
century of direct and indirect cosmic ray detections, the data now spans^l2 
orders of magnitude in energy and about 32 in flux as shown in Figure iB A 
strikingly simple power law spectrum over these many decades in energy is 
observed: up to energies ~ 10^^ eV the spectrum is J{E) oc E ' with 7 ~ 2.7 
and becomes steeper for higher energies, 7 ~ 3. Most of these cosmic rays are 
believed to be accelerated in Galactic Supernovae shocks via Fermi acceleration 
and to propagate diffusively in Galactic magnetic fields. This hypothesis is yet 
to be confirmed from direct observations of Supernova remnants. In addition, 
the energetic requirements at the highest energies call for alternative sources 
which are expected to originate from extragalactic sources since proton Larmor 
radii become larger than Galactic dimensions at these extreme energies. 

Once extragalactic sources of protons are considered, their spectrum should 
not continue as a simple power law due to the effect of the cosmic background 
radiation along their path. In the rest frame of 10^*^ eV protons the cosmic 
microwave background corresponds to gamma rays with energies above the A 
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Figure 1: Cosmic ray spectrum compiled by S. Swordy [2]. 



resonance. Once protons reach these energies, they lose energy very efficiently 
through the production of pions. As ultra-high energy protons travel across 
intergalactic space, the pion photo-production process off the cosmic microwave 
backgrouad quickly reduces their energy to below ~ 7 x 10^^ eV as seen in 
Figure 20 If protons of 10^° eV are observed, their sources are limited to 
distances of less than about 50 Mpc from Earth. This threshold energy loss 
process should introduce a clear feature and a cutoff around ^ iq^Q-^o gy 
the spectrum of extragalactic protons even if the injection spectrum at the 
source is a power law. This photopion production gives rise to the well-known 
Greisen-Zatsepin-KuzmiiH (GZK) cutoff which is shown in Figure 3 by the 
dashed line for a uniform distribution of sources. 

Figure M shows the expected spectrum (flux multiplied by versus en- 
ergy) for extragalactic proton sources (dashed line) and the flux observed by 
the Akeno Giant Air Shower Array ( AGASA) . AGAS A has accumulated many 
hundreds of events (728) with energies above 10^^ eV and 8 events above 10^° 
eVcl As can be seen from the figure, the data shows no indication of a cutoff. 
This surprise has triggered considerable interest and the proposal of a number 
of exotic scenarios designed to explain these data. In addition to the ultra- 
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Figure 2: Proton energy vs. propagation distance through the 2. 7-K cosmic background 
radiation for the indicated initial energies from [2]. 



high energy events detected by AGASA, other experiments such as Fly's Eye, 
Haverah Park, Yakutsk, and Volcano Ranch, have reported a number of events 
at the highest energies (see, e.g., [6] for reviews). 

These observations are surprising because not only the propagation of par- 
ticles at these energies is prone to large energy losses, but the energy require- 
ments for astrophysical sources to accelerate particles to > 10^" eV are extraor- 
dinary. Different primaries do not ease the difficulty. Ultra-high energy nuclei 
are photodisintegrated on shorter distances due to the infrared background 
while the radio background constrains photons to originate from even closer 
systems B 

The shape of the GZK cutoff depends on the source input spectrum and the 
distribution of sources in space as well as in the intergalactic magnetic field. 
In Figure 4, we contrast the flux observed by AGASA with a monte-carlo 
of the expected flux for proton sources distributed homogeneously (shaded 
region) or distributed likfi galaxies (hatched region) with injection spectrum 
J{E) (X and 7 = SEI We modeled the distribution of ultra-high energy 
cosmic ray (UHECR) sources by using the galaxy distribution measured by 
the recent IRAS redshift survey know as PSCzEHl As can be seen from the 
figure, even allowing for the local overdensity the observations are consistently 
above the theoretical expectation. In fact, when we normalize our simulations 
by requiring that the number of events with E > 10^^ eV equals the AGASA 
observations (728), we find that the number of expected events for E > 10^^ 
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Figure 3: E^J{E) versus E observed by AGASA and expected flux for a uniform extragalac- 
tic source distribution [4]- 

eV is only 1.2 ± 1.0 for the PSCz case, i.e., "6 tr" away from the observed 8 
events. 

The gap between observed flux and model predictions narrows as the in- 
jection spectrum of the UHECJ^ sources becomes much harder than 7 = 3. 
For 7 = 2.1 (shown in Figure aa), the number of expected events above 10^° 
eV reaches 3.3 ± 1.6 for a homogeneous distribution while for the PSCz catalog 
it is 3.7 ± 2.0. This trend can be seen also in Figure 6, where mean fluxes for 
7 — 1.5, 2.1 and 2.7 are shown. 

In addition to the presence of events past the GZK cutoff, there has been 
no clear counterparts identified in the arrival direction of the highest energy 
events. If these events are protons or photons, these observations should be 
astronomical, i.e., their arrival directions should be the angular position of 
sources. At these high energies the Galactic and extragalactic magnetic fields 
should not affect proton orbits significantly so that even protons would point 
back to their sources within a few degrees. Protons at 10^" eV propagate 
mainly in straight lines as they traverse the Galaxy since their gyroradii are ~ 
100 kpc in fiG fields whichjs-typical in the Galactic disk. Extragalactic fields 
are expected to be ^ /^Gli^'E3 and induce at most ~ 1° deviation from the 
source. Even if the Local Supercluster has relatively strong fields, the high- 
est energy events may deviate at most ^ 10°E3o At present, no correlations 
between arrival directions and plausible optical counterparts such as sources 
in the Galactic plane, the Local Group, or the Local Supercluster have been 
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Figure 4: Simulated fluxes for the AGASA statistics of 728 events above 10^^ eV, and 7 = 3, 
using a homogeneous source distribution (\ hatches) and the PSCz distribution (dense / 
hatches). The solid and dashed lines are the results of the analytical calculations for the 
same two cases. The dash-dotted and dash- dot- dot- dotted lines trace the mean simulated 
fluxes for the homogeneous and the PSCz cases, (see [9]). 

clearly identified. Ultra high energy cosmic ray data are consistent with an 
isotropic distribution of sources in sharp contrast to the anisotropic distribu- 
tion of light within 50 Mpc from Earth. 

The absence of a GZK cutoff and the isotropy of arrival directions are two 
of the many challenges that models for the origin of UHECRs face. This is an 
exciting open field, with many scenarios being proposed but no clear resolu- 
tion. Not only the origin of these particles may be due to physics beyond the 
standard model of particle physics, but their existence can be used to consteain 
extensions of the standard model such as violations of Lorentz invarianceuS 

In the next section, we discuss the issues involved in the propagation of 
UHECRs from source to Earth. We then summarize the astrophysical Zevatron 
proposals followed by the models that involve new physics. To conclude, future 
observational tests of UHECR models and their implications are discussed. 
Fore recent reviews see [16-19]. 

2 Propagation - Losses and Magnetic Fields 

Before contrasting plausible candidates for UHECR sources with the observed 
spectrum and arrival direction distribution, we discuss the effects of propa- 
gation from source to Earth. Propagation studies involve both the study of 
losses along the primaries' path, such as the photopion production responsible 
for the GZK cutoff, as well as the structure and magnitude of cosmic magnetic 
fields that determine the trajectories of charged primaries and influence the 



Figure 5: Simulated fluxes for the AGASA statistics of 728 events above 10^^ eV, and 
7 = 2.1, using a homogeneous source distribution with Zmax = 0.1 (/ hatches), the PSCz 
distribution with Zmax =0.1 (horizontal hatches), and a homogeneous source distribution 
with Zmax = 1 f\ hatches) (see [9]). 

development of the associated electromagnetic cascadeS 

For primary protons the main loss prpresses are photopion jproduction off 
the CMB that gives rise to the GZK cutofcl and pair productionEil For straight 
line propagation, the cutoff should be present at ^ 5 x 10^^ eV and a significant 
number of hard sources should be located within ^ 50 Mpc from us. As we 
discussed above, even with the small number of accumulated events at the 
highest energies, the AGASA spectrum is incompatible with a GZK cutoff for 
a homogeneous extragalactic source distribution. The shape flt-ttcj cutoff can 
be modified if the distribution of sources is not homogeneoumOcJ and if the 
particle trajectcpes_axfi not rectilinear (e.g., the case of sizeable intergalactic 
magnetic fields) bl-N-y 

Charged particles of energies up to 10^° eV can be deflected significantly 
in cosmic magnetic fields. In a constant magnetic field of strength B — BqiiG, 
particles of energy E — _E2olO^°eV and charge Ze have Larmor radii of — 
110 kpc {E2o/BqZ). If the UHECR primaries are protons, only large scale 
intergalactic magnetic fields affect, their propagation significantly unless the 
Galactic halo has extended fieldsEj For higher pZl,_tJie Galactic magnetic field 
can strongly affect the trajectories of primariesEj'EJ 

Whereas Galactic magnetic fielda,are reasonably well studied, extragalactic 
fields are still very ill understoodLj Faraday rotation measures show large 
magnitude fields ^G) in the central regions of clusters of galaxies. In regions 
between clusters, the presence of magnetic fields is evidenced by synchrotron 
emission but the strength and structure are yet to be determined. On the 
largest scales, limits can be imposed by the observed isotropy of the CMB 
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Figure 6: Propagated spectrum for source spectral index of -y 



= 1.5, 2.1, 2.7. 



and by a statistical interpretation of Faraday rotation measures of light from 
distant quasars. The isotrop^ppf the CMB can constrain the present horizon 
scale fields B^j-i^Z x 10"^ Gt3 Although the distribution of Faraday rotation 
measures have large non-gaussian tails, a reasonable limit can be derived using 
the median of the distribution in an inhomogcneous universe: for fields assumed 
to be constant on the present horizon scale, i3„-i< 10^^ G; for fields with 50 

Mpc coherencCjIength, SsoMpc^S x 10~^ G; while for 1 Mpc coherence length, 
Bmpc^ 10~* GI13 These limits apply to a Qbh'^ = 0.02 universe and use quasars 
up to redshift z — 2.5. Local structures can have fields above these upper limits 



as . 
2.5 



Lg^as they are not common along random lines of site between z — and 

13113 



Of particular interest is the field in the local 10 Mpc volume around us. If 
the Local Supercluster has fields of about 10^^ G or larger, the propagation of 
ultra high energy protons could become diffusive and the spectrum ancLapgiib 
lar distribution at the highest energies would be significantly modifiedtj'ESEj 
For example, in Figure 7E3 a source with spectral index 7^2 that can reach 
Emax^^Q^'^ eV is constrained by the overproduction of lower energy events 
around 1 to 10 EeV (EeV = 10^* eV). Furthepnore, the structure and mag- 
nitude of magnetic fields in the Galactic halcj'Ej or in a possible Galactic 
wind can also affect the observed UHECRs. In particular, if our Galaxy has 
a strong magnetized wind, what appears to be an isotropic distribution in ar- 
rival directions, may have originated on a small region of the sky such as the 
Virgo clustered In the future, as sources oLUHECRs are identified, large scale 
magnetic fields will be better constrainedcil 

If cosmic rays are heavier nuclei, the attenuation length is shorter than that 
for protons due to photodisintegration on the infrared backgroundQ However, 
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Figure 7: Flux vs. Energy with Emax = 10^^ eV at source. Choices of source distance 
r(Mpc), spectral index 7, proton luminosity Lp(erg/s), and LSC field B(/iG) are: solid line 
(13 Mpc, 2.1, 2.2 X 10*3 erg/s, 0.05 fiG); dotted line (10 Mpc, 2.1, lO^^ erg/s, 0.1 fiG); 
dashed Une (10 Mpc, 2.4, 3.2 xlQ-^^ erg/s, 0.1 fiG); and dashed-dotted Une (17 Mpc, 2.1, 
3.3 X 10^^ erg/s, 0.05 fiG). 

UHE nuclei may be of Galactic origin. For large enough charge, the trajeiitories 
of UHE nuclei are significantly affected by_the Galactic magnetic fielcEJ such 
that a Galactic origin can appear isotropics The magnetically induced distor- 
tion of the flux map of UHE events can give rise to some higher flux-tcgions 
where caustics form and some much lower flux regions (blind spots) Such 
propagation effects are one of the reasons why full-sky coverage is necessary 
for resolving the UHECR puzzle. 

The trajectories of photon primaries are not affected by magnetic fields, 
but energy losses due to tJae radio background constrain photons to originate 
from systems at <10MpcEj If associated with luminous systems, sources of 
UHE photons should point back to their nearby sources. The lack of coun- 
terpart identifications suggests that if the primaries are photons, their origin 
involves physics beyond the standard model. 

3 Astrophysical Zevatrons 

The puzzle presented by the observations of cosmic rays above 10^° eV have 
generated a number of proposals that wc divide here as Astrophysical Zeva- 
trons and New Physics models. Astrophysical Zevatrons are also referred to 
as bottom-up models and involve searching for acceleration sites in known as- 
trophysical objects that can reach ZeV energies. New Physics proposals can 
be either hybrid or pure top-down models. Hybrid models involve Zevatrons 
and extensions of the particle physics standard model while top-down models 
involve the decay of very high mass relics from the early universe and physics 



way beyond the standard model. Here we discuss astrophysical Zevatrons while 
new physics models are discussed in the next section. 

Cosmic rays can be accelerated in astrophysical plasmas when large-scale 
macroscopic motions, such as shocks, winds, and turbulent flows, are trans- 
ferred to individual particles. The maximum energy of accelerated particles, 
Emax, can be estimated by requiring that the gyroradius of the particle be 
contained in the acceleration region: i?max = ZeB L, where Ze is the charge 
of the particle, B is the strength and L the coherence length of the magnetic 
field embedded in the plasma. For fi'max^lO^'^ eV and Z ^ 1, the only known 
astrophysical sources with reasonable BL products are neutron stars {B ^ 10^^ 
G, L - 10 km), active galactic nuclei (AGNs) (B ~ 10^ G, L ~ 10 AU), radio 
lobes of AGNs {B ~ O.l/xG, L ~ 10 kpc), and clusters of galaxies {B ~ /xG, 
L ^ 100 kpc). In Figure 8, we highlight the B vs. L for objects that can reach 
Emax = 10^° eV with Z = 1 (dashed hne) and Z = 26 (sohd hne). We discuss 
each of these candidates below. 
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Figure 8: B vs. L, for Emax = 10^^ eV, Z = 1 (dashed line) and Z = 26 (solid line) from 
[16]. 

Clusters of Galaxies: Cluster shocks are reasonable sites to consider for 
UHECR acceleration, since particles with energy up to -Emax can be contained 
by cluster fields. However, efficient losses due to photopion production off the 
CMB during the propagaticm^inside the cluster limit UHECRs in cluster shocks 
to reach at most 10 EcVlj 

AGN Radio Lofees/JNFext on the list of plausible Zevatrons are extremely 
powerful radio galaxies eZI Jets from the central black- hole of an active galaxy 
end at a termination shock where the interaction of the jet with the intergalac- 
tic medium forms radio lobes and 'hot spots'. Of special interest are the most 
powerful AGNs where shocks can accelerate particles to energies well above an 
EeV via the first-order Fermi mechanism. These sources may be responsible 
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for the flux of UHECRs up to the GZK cutoflS 

A nearby specially powerful source may be able to reach energies past the 
cutoff. However, extremely powerful AGNs with radio lobes and hot spots are 
rare and far apart. The closest known object is M87 in the Virgo cluster 
18 Mpc away) and could be a main source of UHECRs. Although a single 
nearby source with especially hard spectra may fit the spectrum for a given 
strength and structure of the intergalactic magnetic fieldpa it is unlikely to 
match the observed arrival direction distribution. If M87 is the primary source 
of UHECRs a concentration of events in the direction of M87 or the Virgo 
cluster should be seen in arrival directions. No such hot spot is observed {Hot 
spot in Table 1). The next known nearby source after M87 is NGC315 which 
is already too far at a distance of ~ 80 Mpc. Any unknown source between 
M87 and NGC315 would likely contribute a second hot spot, not an isotropic 
distribution. The very distant radio lobes will contribute a GZK cut spectrum 
which is not observed. 

The lack of a clear hot spot in the direction of M87 has encouraged the 
idea that a strong Galactic magnetic wind may exist that could help isotropize 
the arrival directions q£ UHECRs. A Galactic wind with a strongly magnetized 
azimuthal componenio {GalacticWind in Table 1) can significantly alter the 
paths of UHECRs such that the observed arrival directions of events above 10^° 
eV would trace back to the North Galactic Pole which is close to the Virgo 
where M87 resides. If our Galaxy has such a wind is yet to be determined. The 
proposed wind would focus most observed events into the northern Galactic 
pole and render point source identification fruitlesstj Future observations of 
UHECRs from the Southern Hemisphere by the Southern Auger Observatory 
will provide precious data on previously unobserved parts of the sky and help 
distinguish plausible proposals for the effect of local magnetic fields on arrival 
directions. Full sky coverage is a key discriminator of such proposals. 

AGN - Central Regions: The powerful engines that give rise to the observed 
jets and radio lobes are located in the central regions of active galaxies and are 
powered by the accretion of matter onto supermassive black holes. It is rea^ 
sonable to consider the central engines themselves as the likely acceleratorsc2l 
In principle, the nuclei of generic active galaxies (not only the ones with radio 
lobes) can accelerate particles via a unipolar inductor not unlike the one oper- 
ating in pulsars. In the case of AGNs, the magnetic field may be provided by 
the infalling matter and the spinning black hole horizon provides the imperfect 
conductor for the unipolar induction. 

The problem with AGNs as UHECR sources is two-fold: first, UHE par- 
ticles face debilitating losses in the acceleration region due to the intense ra- 
diation field present in AGNs, and second, the spatial distribution of objects 
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should give rise to a GZK cutoff of the observed spectrum. In the central 
regions of AGNs, loss processes are expected to downgrade particle energies 
well below the maximum achievable energy. This limitation has led to the 
proposal that quasar remnants, supermassive black. holes in centers of inactive 
galaxies, are more effective UHECR accelerators EJ In this case, losses at the 
source are not as significant but the propagation from source to us should still 
lead to a clear GZK cutoff since sources would be associated with the large 
scale structure of the galaxy distribution {LSS in Table 1). From Figure 4-6, 
these models can only succeed if the source spectrum is fairly hard (7<2)lI 
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Neutron Stars: Another astrophysical system capable of accelerating UHE- 
CRs is a neutron star. In addition to having the ability to confine lO^'^eV 
protons (Figure 8), the rotation energy of young^ neutron stars is more than 
sufficient to match the observed UHECR fluxesEJ However, acceleration pro- 
cesses inside the neutron star light cylinder are bound to fail much like the 
AGN central region case: ambient magnetic and radiation fields induce signifi- 
cant losses. However, the plasma that expands beyond the light cylinder is free 
from the main loss processes and may be accelerated to ultra high energies. 

One possible source of UHECR past the GZK cutoff is the early evolution of 
neutron stars. In particular, newly formed, rapidly rotating neutron stars may 
accelerate iron-nuclei to UHEs through relativistic MHD winds beyond their 
light cylinderscS This mechanism naturally leads to very hard injection spectra 
(7 ~ 1) (see Table 1). As seen in Figure 6, 7 1 improves the agreement 
between predicted flux and observations for energies above 10^" eV. In this 
case, UHECRs originate mostly in the Galaxy and the arrival directions require 
that the primaries be heavier nuclei. Depending on the structure of Galactic 
magnetic fields, the trajectories of iron nuclei from Galactic neutron stars ma^ 
be consistent with the observed arrival directions of the highest energy eventstj 
Moreover, if cosmic rays of a few times 10^* eV are protons of Galactic origin, 
the isotropic distribution observed at these energies is indicative of the diffusive 
effect of the Galactic magnetic fields on iron at ~ 10^° eV. This proposal should 
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be constrained once the primary composition is clearly determined (see Iron 
in Table 1). 

It has also been suggested that young extragalactic highly magnetized neu- 
tron stars (magnetars) onay be sources of UHE protons which are accelerated 
by reconnection eventscj These would be prone to a GZK cut spectrum and 
would need a very hard injection spectrum to become viable explanations. 

Gamma-Ray Bursts: Transient high energy phenomena such as gaipipa- 
ray bursts (GRBs) may also be a source of ultra-high energies protonsc3 In 
addition to both phenomena having unknown origins, GRBs and UHECRs have 
other similarities that may argue for a common source. Like UHECRs, GRBs 
are distributed isotropically in the sky, and the average rate of 7-ray energy 
emitted by GRBs is comparable to the energy generation rate of UHECRs 
of energy > 10^^ eV in a redshift independent cosmological distribution of 
sources, both have « 10^'*erg /Mpc^/yr. 

However, recent GRB counterpart identifications argue for a strong cos- 
mological evolution for GRBs. The redshift dependence of GRB distribution 
is such that the flux of UHECR ass£iciated with nearby GRBs would be too 
small to fit the UHECR observationsEJ In addition, the distribution of UHECR 
arrival directions and arrival times argues against the GRB-UHECR common 
origin. Events past the GZK cutoff require that only GRBs from <50 Mpc con- 
tribute. Since less than about one burst is expected to have occurred within 
this region over a period of 100 yr, the unique source would appear as a con- 
centration of UHECR events in a small part of the sky (a Hot spot in Table 
1). In addition, the signal would be very narrow in energy /S.E/ E ~ 1. Again, 
a strong intergalactic magnetic field can ease some of these difficulties giving 
a very large dispersion in the arrival tinie, and direction of protons produced 
in a single burst {large Bjgm in Table l)ua Finally, if the observed small scale 
clustering of arrival directions is confirmed by future experiments with clus- 
ters having some lowet-energy events clearly precede higher energy ones, bursts 
would be invalidatedCJ 

4 New Physics Models 

The UHECR puzzle has inspired a number of different models that involve 
physics beyond the standard model of particle physics. New Physics proposals 
can be top-down models or a hybrid of astrophysical Zevatrons with new par- 
ticles. Top-down models involve the decay of very high mass relics that could 
have been formed in the early universe. 

The most economical among hybrid proposals involves a familiar exten- 
sion of the standard model, namely, neutrino masses. If some flavor of neu- 
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trinos have mass (e.g., ^ O.leV), the rehc neutrino background is a target for 
extremely high energy neutrinos to interact, and generate other particles by 
forming a Z-boson that subsequently decayEll (see v Z burst in Table 2). If the 
universe has very luminous sources (Zevatrons) of extremely high energy neu- 
trinos (^ 10^^ eV), these neutrinos would traverse very large distances before 
annihilating with neutrinos in the smooth cosmic neutrino background. UHE 
neutrino Zevatrons can be much further than the GZK limited volume, since 
neutrinos do not suffer the GZK losses. But if the interaction occurs through- 
out a large volume, the GZK feature should also be observed. For plausible 
neutrino masses ~ O.leV, the neutrino background is very unclustered, so the 
arrival direction for events should be isotropic and small scale clustering may 
be a strong challenge for this proposal. The weakest link in this proposal is 
the nature of a Zevatron powerful enough to accelerate protons above tens of 
ZeVs that can produce neutrinos as secondaries with ^ ZeV. This Zevatron is 
quite spectacular, requiring an energy generation in excess of presently known 
highest energy sources (referred to as ZeV v's in Table 2). 

Another suggestion is that the UHECR primary is a new hadronic particle 
that is also accelerated in Zevatrons. The mass of a hypothetical hadronic 
primary can be limited by the shower development of the Fly's Eye highest 
energy event to be below ^50 GeVeS As in the Z-burst proposal, a neutral 
particle is usually harder to accelerate and are usually created as secondaries 
of even higher energy charged primaries. But once formed these can traverse 
large distances without being affected by cosmic magnetic fields. Thus, a 
signature for future experiments of hybrid models that invoke new particles 
as primaries is a clear correlation between the position of powerful Zevatrons 
in the sky sudi as distant compact radio quasars and the arrival direction of 
UHE eveatsca Preliminary evidence for such a correlation has been recently 
report edE3 

Another exotic primary that can be accelerated to ultra high energies by 
astrophysical systems is the vorton. Vortons are small loops of superconducting 
cosmic string stabilized by the angular momentum of charge carriersEiJ Vortons 
can be a component of the darkjuatter in galactic halos and be accelerated 
by astrophysical magnetic fieldsEj Vortons as primaries can be constrained by 
the observed shower development profile. 

It is possible that none of the astrophysical scenarios or the hybrid new 
physics models are able to explain present and future UHECR data. In that 
case, the alternative is to consider top-down models. Top-down models involve 
the decay of monopole-antimonoploe pairsp ordinary and superconducting 
cosmic strings, cosmic necklaces, and superheavy long-lived relic particles. The 
idea behind these models is that relics of the very early universe, topological 
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defects (TDs) or superheavy relic (SHR) particles, produced after or at the end 
of inflation, can decay today and generate UHECRs. Defects, such as cosmic 
strings, domain walls, and magnetic monopoles, can be generated through the 
Kibble mechanism as symmetries are broken with the expansion and cooling 
of the universe. Topologically stable defects can survive to the present and 
decompose into their constituent fields as they collapse, annihilate, or reach 
critical current in the case of superconducting cosmic strings. The decay prod- 
ucts, superheavy gauge and higgs bosons, decay into jets of hadrons, mostly 
pions. Pions in the jets subsequently decay into 7-rays, electrons, and neutri- 
nos. Only a few percent of the hadrons are expected to be nucleons. Typical 
features of these scenarios are a predominant release of 7-rays and neutrinos 
and a QCD fragmentation spectrum which is considerably harder than the case 
of Zevatron shock acceleration. 

ZeV energies are not a challenge for top-down models since symmetry 
breaking scales at the end of inflation typically are ^ 10^"'^ eV (typical X- 
particle masses vary between j^q^^"^^ eV). Fitting the observed flux of 
UHECRs is the real challenge since the typical distances between TDs is the 
Horizon scale, i. e., about several Gpc. The low flux hurts proposals based on 
ordinary and superconducting cosmic strings which are distributed throughout 
space [Extragal. TD in Table 2). Monopoles usually suffer the opposite prob- 
lem, they would in general be too numerous. Inflation succeeds in diluting the 
number density of monopoles and makes them too rare for UHECR produc- 
tion. To reach the observed UHECR flux, monopole models usually involve 
some degree of fine tuning. If enough monopoles and antimonopoles survive 
from the early universe, they may form a bound state, named monopolonium, 
that can decay generating UHECRs. The lifetime of monopolonia majj_be too 
short for this scenario to succeed unless they are connected by strings Efl 

Once two symmetry breaking scales are invoked, a combination of hori- 
zon scales gives room to reasonable number densities. This can be arranged 
for cosmic strings that end in monopoles jnaking a monopole string network 
or even more clearly for cosmic necklaces E3 Cosmic necklaces are hybrid de- 
fects where each monopole is connected to two strings resembling beads on a 
cosmic string necklace. Necklace networks may evolve to configurations that 
can fit the UHECR flux which is ultimately generatedLhji. the annihilation of 
monopoles with antimonopoles trapped in the stringE3u In these scenarios, 
protons dominate the flux in the lower energy side of the GZK cutoff while 
photons tend to dominate at higher energies depending on the radio back- 
ground (see Figure 9 and Extragal. TD in Table 2). If future data can settle 
the composition of UHECRs from 0.01 to 1 ZeV, these models can be well 
constrained. In addition to fitting the UHECR flux, topological defect models 
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Figure 9; Proton and 'y-ray fluxes from necklaces for mx = 10^'* GeV (dashed lines), 10^^ 
GeV (dotted lines), and 10^® GeV (solid lines) normalized to the observed data, 'y-high and 
■y-low correspond to two extreme cases of ^ -ray absorption (see, [56]). 



are constrained by limits on the flux of high energy photons, from 10 MeV to 
100 GeV, observed by EGRET. 




loo[E(„V)] 



Figure 10: SHRs or monopolia decay fluxes (for mx = 10^* GeV): nucleons from the 
halo (protons), -y -rays from the halo (gammas) and extragalactic protons. Solid, dotted and 
dashed curves correspond to different model parameters (see [56]). 

Another interesting possibility is the recent proposal that UHECRs are 
produced by the decay of unstable superheavy relics that live much longer 
than the age of the universeEll SHRs may be produced at the end of infla- 
tion by non-thermal effects such as a varying gravitational field, parametric 
resonances during preheating, instant preheating, or the decay of topological 
defects. These models need to invoke special symmetries to insure unusually 
long lifetimes for SHRs and that a sufficiently small percentage decays today 
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producing UHECRsEZlH As in the topological defects case, the decay of these 
relics also generates jets of hadrons. These particles behave like cold dark mat- 
ter and could constitute a fair fraction of the halo of our Galaxy. Therefore, 
their halo decay products would not be limited by the GZK cutoff allowing 
for a large flux at UHEs (see Figure 10 and SHRs in Table 2). Similar sig- 
natures can occur if topological defects are microscopic, such as monopolonia 
and vortons, and decay in the Halo of our Galaxy {Halo TD in Table 2). In 
both cases [SHRs and Halo TD) the composition of the primary would be a 
good discriminant since the decay products are usually dominated by photons. 



Table 2: NEW PHYSICS 

v Z burst Extragal TD 



Halo TD 



SHRs 



Composition 
Source 7 
Sky Distrib. 
Problem:Pho- 



Photon 
Z frag 
Isotropic 
ZeV v's 



Pho+GZK p 
QCD frag 
Isotropic 
flux 



Photon 
QCD frag 
Gal Halo 
origin 



Photon 
QCD frag 
Gal Halo 
lifetime 



Future experiments should be able to probe these hypotheses. For in- 
stance, in the case of SHR and monopolonium decays, the arrival direction 
distribution should be close to isotropic but show an asymmetry due to the 
position of the Earth in the Galactic HalcE^ and the clustering due to small 
scale dark matter inhomogeneitiesE^I Studying plausible halo models for their 
expected asymmetry and inhomogeneitis will help constrain halo distributions 
especially when larger data sets are available in the future. High energy gamma 
ray experiments such as GLASTLwill also help constrain SHR models via the 
electromagnetic decay productsEd 



5 Conclusion 

Next generation experiments such as the High Resolution Fly's Eye which 
recently started operating, the Pierre Auger Project which is now under con- 
struction, the proposed Telescope Array, and the EUSO and OWL satellites 
will significantly improve the data at the extremely-high end of the cosmic ray 
spectrumB With these observatories a clear determination of the spectrum and 
spatial distribution of UHECR sources is within reach. 

The lack of a GZK cutoff should become clear with HiRes and Auger and 
most extragalactic Zevatrons may be ruled out. The observed spectrum will 
distinguish Zevatrons from new physics models by testing the hardness of the 
spectrum and the effect of propagation. Figure 11 shows how clearly Auger 
will test the spectrum in spite of their clustering properties. The cosmography 
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Figure 11: Simulated fluxes for the Auger projected statistics of 9075 events above 10^® eV, 
and 7 = 3, using a homogeneous source distribution (\ hatches) and the PSCz distribution 
(/ hatches). The solid and dashed lines are the results of the analytical calculations for the 
same two cases. The dash-dotted and dash- dot- dot- dotted lines trace the mean simulated 
fluxes for the homogeneous and the PSCz cases, (see [7]). 



of sources should also become clear and able to discriminate between plausible 
populations for UHECR sources. The correlation of arrival directions for events 
with energies above 10^° eV with some known structure such as the Galaxy, 
the Galactic halo, the Local Group or the Local Supercluster would be key in 
differentiating between different models. For instance, a correlation with the 
Galactic center and disk should become apparent at extremely high energies 
for the case of young neutron star winds, while a correlation with the large 
scale galaxy distribution should become clear for the case of quasar remnants. 
If SHRs or monopolonia are responsible for UHECR production, the arrival 
directions should correlate with the dark matter distribution and show the 
halo asymmetry. For these signatures to be tested, full sky coverage is essen- 
tial. Finally, an excellent discriminator would be an unambiguous composition 
determination of the primaries. In general. Galactic disk models invoke iron 
nuclei to be consistent with the isotropic distribution, extragalactic Zevatrons 
tend to favor proton primaries, while photon primaries are more common for 
early universe relics. The hybrid detector of the Auger Project should help 
determine the composition by measuring the depth of shower maximum and 
the muon content of the same shower. The prospect of testing extremely high 
energy physics as well as solving the UHECR mystery awaits improved ob- 
servations that should be coming in the next decade with experiments under 
constmiition such as AugefiiLpr in the planning stages such as the Telescope 
ArrayP EUSOP and OWlN 
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